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ABSTRACT 

We analyse the kinematics of disc stars observed by the RAVE survey in and beyond the Solar 
neighbourhood. We detect significant overdensities in the velocity distributions using a tech- 
nique based on the wavelet transform. We find that the main local kinematic groups are large 
scale features, surviving at least up to ~ 1 kpc from the Sun in the direction of anti-rotation, 
and also at ~ 700 pc below the Galactic plane. We also find that for regions located at different 
radii than the Sun, the known groups appear shifted in the vr-v$ velocity plane. For example, 
the Hercules group has a larger azimuthal velocity for regions inside the Solar circle and a 
lower value outside. We have also discovered a new group at (U, V) — (92, —22) kins" 1 
in the Solar neighbourhood and confirmed the significance of other previously found groups. 
Some of these trends detected for the first time are consistent with dynamical models of the 
effects of the bar and the spiral arms. More modelling is required to definitively characterise 
the non-axisymmetric components of our Galaxy using these groups. 

Key words: Galaxy: kinematics and dynamics - Galaxy: disc solar neighbourhood - Galaxy: 
structure - Galaxy: evolution - 



1 INTRODUCTION 



The presence of kinematic structures in the velocity distribution 
of the Solar neighbourhood has been well esta blished through 
the analysis of H i pparcos and other surveys dDehnenl Il998t 
iFamaev et ai]|2005r . lAntoia et alj|2008l . hereafter A08). Some of 
the structures may be related to orbital effects of the bar dDehnenl 
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2011L hereafter All) or both 
1 120 111 All). Thus, they can pro- 
vide constraints on the properties of the non-axisymmetries of the 
Milky Way (MW) disc. However, the origin of each structure is not 
yet clear. Some of them could be remnants of disrupted disc stellar 
clusters or could b e associated to accretion events (see review in 
lAntoia et al1l2010h - 



The RAVE (RAdial Velocity Experiment, Steinm etz et alj 

2006) spectroscopic survey provides the means for a giant leap for- 
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Figure 1. Spatial distribution in X-Y and X-Z of RAVE data with 6D 
phase space coordinates (black points). Grey points are stars with \Z\ ^ 
300 pc. Orange points are stars in the GCSIII. Other colours are the RAVE 
subsamples defined in the text. The black dashed line shows the Solar circle. 



Table 1. Properties of the subsamples defined in the text. For each subsam- 
ple we give the number of stars, median distance in kpc, median relative 
error in distance, and median errors in heliocentric velocities (for local sub- 
samples) or in cylindrical velocities (for distant subsamples ) in km s . 
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ward in the understanding of the origin of the kinematic groups, 
allowing us to answer primary questions such as: what are the sig- 
nificant groups in the Solar volume? are they local features? do they 
change when we move far from the Sun? The RAVE Solar neigh- 
bourhood sample contains already more stars than any previously 
available k inematic data set, such as the Geneva-Copenhagen sur- 
vey (GCS. lHolmberg et alj ]2009). But more importantly, the vol- 
ume sampled by RAVE is much larger, giving kinematics represen- 
tative of distant regions of the thin and thick discs. In this Letter, 
we use the RAVE data to revisit the kinematic groups in the Solar 
neighbourhood and trace them, for the first time, far from the Sun. 



2 THE RAVE DATA AND METHODS 

RAVE is a magnitude limited multi-fibre spectroscopic survey in 
the range 9 < / < 13 using the 1.2-m UK Schmidt Telescope 
of the Angl o-Australian Obser vatory (AAO). We use here internal 
data from Burne tt et al. I J20T l]), w ho derived distances for stars in 
the sec ond and thir d data releases dZ witter et al.l2008llSiebert et alj 
l2011bl) . Following lBurnettetalJ fa) 111) , we increase the distances 



by 10% for hot dwarfs and decrease them by the same factor for gi- 
ants. The data set i ncludes also pr oper motions mainly f rom PPMX 
dRoser et alj|2008l) and UCAC2 dZacharias et"al]|2004 . making a 
sample of 202843 stars with 6D phase space information. 

The spatial distribution of the sample is shown in Fig.Q](black 
points). We use the Galactic Cartesian coordinate system with the 
Su n at X = —7.8 kpc, Y = 0. kpc (second mod el in table 2 
of iMcMillan & Binned 120101) and Z = 0.014 kpc dBinnev et al] 
1 19971) . The general spatial distribution can be explained by the fact 
that RAVE is a southern hemisphere survey (S < 2°). The cone 
through negative Y is due to the additional colour cut J — K > 0.5 
that was applied for latitudes \b\ < 25° and longitudes 230° < 
I < 315° to select giant stars and probe larger distances in the 
disc. There are 112808 stars with \Z\ ^ 300 pc (grey points). As 
a comparison, throughout the paper we will us e the G CS sample 
(hereafter, GCSIII) taken from iHolmberg et al.l d2009l) . The main 
differences between the two data sets are that RAVE is a southern 
hemisphere survey whereas GCSIII covers the whole sky, that it is 
much fainter and that it is not restricted to F and G dwarfs. We plot 
the GCSIII with orange points in Fig. [T] The volume sampled by 
the RAVE stars is much larger and the kinematics of the thin disc 
can be now traced beyond ~ 1 kpc in the direction of anti-rotation. 

We study the kinematics of stars in several spatially distinct 
samples whose main properties are in Table Q] First, we define the 
local subsample (LS, blue points in Fig.QJ as those stars with \Z\ ^ 
300 pc and in-plane distances from the Sun dcos(fo) ^ 200 pc. 
There are 57178 stars in the LS, which is more than four times the 
number in the GCSIII. Secondly, we select stars in three regions 
in the plane (\Z\ ^ 300 pc) corresponding to cubic volumes of 
0.6 kpc size at a distance of ~ 1 kpc from the Sun in the range 

— 1.3 ^ Y < — 0.7kpc but at different Galactocentric radius: 
i) Sr : Solar radius —8.1 < X ^ —7.5 kpc (green points), ii) 
Si n : inner radius —7.5 < X ^ —6.9 kpc (red points), and iii) 
Sout: outer radius —8.7 < X —8.1 kpc (purple points). Fi- 
nally, the subsample Z57 (pink points) contains stars in a cylinder 
with dcos(6) ^ 400 pc and height below the plane —700 ^ Z ^ 

— 500 pc. Note that velocity errors in the RAVE subsamples (last 
columns in Table [TJ are dominated by errors in the transverse mo- 
tion (for which distances and proper motions are needed) as the 
radial velocity accuracy is considerably better (et, !os ~ 2 km s _1 ). 

To detect structures in the 2d velocity distributions of these 
subsa mples we use a tec hnique based on the wavelet transform 
(WT) dStarck & Murtaghll2002l ) The WT techniques were first ap- 
plied to detect local str eams by Chereul et al.l dl999h . We use here 
the a trous algorithm dStarck & Murtaghl |2002|) . This consists of 
first applying a set of smoothing filters with progressively increas- 
ing scales to the original data. Then the difference between two suc- 
cessive smoothed planes gives the WT coefficients, which contain 
the information on the structures present in between the two scales 
or sizes. For every position in the velocity plane, a WT coefficient 
equal to corresponds to a constant signal, while a positive (neg- 
ative) value is associated to overdensities (underdensities). Typical 
sizes of the kinematic structures are ~ lOkms" 1 , so in the ab- 
sence of errors they will be prominent on maps in between scales 
of ~ 6 and ~ 15kms . However, in practice the velocity errors 
may determine the scales that will show reliable groups. 

Then we search for maxima among the WT coefficients to de- 
tect the different groups. Fi nally, to anal y se the ir significance we 
use a method developed in ISlezaketalJd 19931) especially suited 
for low intensity levels. This evaluates the probability that the WT 
coefficients are not due to Poisson noise by comparing, in wavelet 
space, the local distribution of stars to what would result for a uni- 
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Figure 2. a) Heliocentric velocities of the local sample (LS). b) Velocity structures in the local sample (LS) obtained with the wavelet transform at scales 
of 6-13 kms -1 . The red solid (dotted) contour shows the 3-cr confidence limit for positive wavelet coefficients or overdensities (negative coefficients or 
underdensities). Red (orange) crosses show the maxima of the detected peaks inside the 3-cr (2-cr) confidence limit, c) Same as b, but for GCSIII sample. 



Table 2. Kinematic groups detected in the local sample (LS) at scales of 6- 
13 kms -1 . The columns give: order in WT coefficient, name of the group 
in previous literature (if any), U velocity of the peak, error in the mean 
U of the stars in the group, V velocity, error in the mean V, approximate 
number of stars and statistical significance. The last three columns show the 
equivalent groups for the GCSIII. Velocities are in km s . 
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form distribution. We will consider coefficients that have probabil- 
ity of being real detections P ^ 99.7% and 95.4% < P < 99.7%, 
similar to the 3-cr and 2-cr significance levels in the Gaussian case, 
respectively. To perform the calculations we use the MR software 
developed by CEA (Saclay, France) and Nice Observatory. An es- 
timate of the number of stars in each group is obtained by counting 
stars enclosed in a circle with the radius of the minimum scale of 
the considered range of scales. 



3 REVISITED SOLAR VOLUME KINEMATICS 

Fig-Ek shows the velocities of stars in the LS. We use a heliocentric 
system U, V, W, with the U towards the Galactic Centre. We focus 
only on the U and V components as th e moving groups, except the 
very young ones l lChereuletal1ll999h . have almost indistinguish- 



able mean W dDehned[l99ll ; ISeabroke & Gilmorell2007l) . Fig.|2p 
corresponds to the positive WT coefficients (overdensities) in the 
LS velocity distribution at scales of 6-13 kms -1 (where 6 and 

13 kms -1 are the scales of the smoothing filters used in the com- 
putation of the WT). This distribution is dominated by clear over- 
densities. Some structures are elongated in V and slightly tilted. We 
find 9 peaks inside the 3-cr confidence level region (marked in red) 
and 10 more that are 2-cr significant (orange). The peaks are num- 
bered from maximum to minimum WT coefficient and are listed in 
Tabled Fig. lit shows the same WT for the GCSIII. We thus find 
similar significant groups in the t wo independent sam ples. But in 
comparison to previous work (e.g. iFamaev et alj |2005. A08), with 
RAVE we have better statistics. 

The four main groups that we find in the LS (from 1 
to 4) are the well-known Coma Berenices, Hyades, Sirius and 
Pleiades. Their positions agree well with the GCSIII (separation 
between maxima < 10 kms -1 ). The only exception is the Hyades 
stream: the stars of the Hyades open cluster in the GCSIII at (U, 
V)=(— 43,— 18) kms -1 (absent in the LS) shift the peak to more 
negative U. The different ranking of these 4 groups in the LS and 
GCSIII may be due to differences in the populations of the samples. 

The Hercules structure is the elongated feature in the U direc- 
tion at V ~ —50 kms -1 . In the LS, group 6 is its main overden- 
sity. Also groups 8, 10 and 13 (this with higher significance than for 
the GCSIII) could be associated to it. Notice that the detection of 
localised maxima in elongated and extended structures such as Her- 
cules may be a bias due to the isotropy of the WT. The red dotted 
contours between the main four groups and the Hercules structure 
correspond to a 3-cr significant underdensity of stars (negative WT 
coefficients ). Another remark able 3-cr underdensity is the so-called 
edge line dSkulianetalJll999h in the upper part of the distribution 
which is present both in the LS and the GCSIII. 

The fifth group in order of WT coefficient (before the H ercules 
stream ) is close to the classic group Wolf 630 studied in Eggen 
(1971). Group 7 is linked to this group in a branch-like shape (here- 
after Wolf elongation). We also find other grou ps that can be associ- 
ated to the ones studied by Eggen] l ll97lL[l99^) . For instance, group 
10 is similar to elndi, group 9 and 1 1 are close to 7Leo, and group 

14 could be identified with 7/Cep. 

Group 12 is a new 3-cr significant group at (U, V)=(92, 
— 22) kms -1 . It is accompanied by group 18 but only at 2-cr level. 
It has around 50 stars but its location in an empty region of the 
distribution makes it conspicuous. Notice that this is a group on 
a highly eccentric orbit, moving towards the Galactic Centre at 
~ 100 kms -1 . There are other groups (15-19) yet to be confirmed 
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Figure 3. Velocity structures at scales of 1 1-22 km s 1 in cylindrical velocities for the subsamples S ou t (outside Solar circle), Sr q (Solar circle), Si n (inside 
Solar circle), and Z57 (—700 ^ Z ^ —500 pc). Yellow squares show the positions of the local sample (LS) groups. Colour code is the same as in Fig.|2p. 



as they are only 2-a significant. Some of them could be tentatively 
related to other known groups. For instance, groups 14, 17 or 1 9 
could be related to Arcturus at V ~ —110 kms' - 1 ( lEggenlll996h , 
We leave a more thorough comparison to gro ups in other studies 
(e.g. lKlement et aljl201lfl Williams et alj|20llh for future work. 



4 STRUCTURES BEYOND THE LOCAL VOLUME 

The three first panels of Fig. [3] show the velocity structures of the 
subsamples Sout, Sr„ and Si n . Now we use a Galactic cylindrical 
coordinate system, vr and v^, which is more suitable for distant 
regions. We use Local Standard of Rest (LSR) circular velocity of 
no = 247 km s -1 and Solar movement with respect to the LSR 
of (U e ,V e ) = (10.0, 11 .0) kms" 1 (second model in Table 2 of 
McMillan & Binnev 2010). For convenience, we take vr and v$ to 
be positive towards the Galactic Centre and towards rotation (clock- 
wise), respectively, to get the same orientation as U—V. Larger WT 
scales (11-22 kms" 1 ) are considered now due to the larger velocity 
errors in these subsamples (Table[TJ. We find that for Sr q and S ou t 
the general shape of the distribution is preserved and some groups 
show a relative configuration similar to the LS. Although we show 
here only three distant regions, the evolution in the location and in- 
tensity of the main peaks as one moves away from the LS to these 
volumes is found to be smooth and continuous. For example, the 
most populated group in Sr (group 1 in the second panel) is the 
Pleiades group. Also groups 2 and 3 are in the Sirius and Hyades 
kinematic positions, respectively. Group 4 is the Hercules stream 
and group 5 corresponds to 7Leo. We also see the Wolf elongation, 
but with no associated maxima. Similar identifications can be done 
for Sout- By contrast, the velocity distribution in Si n looks quite 
different to the LS, its overall shape being elongated with a posi- 
tive slope in the vr direction. There are several central groups but 
they are distributed along the general elongated structure and do 
not show the same relative configuration as for LS, Sr q and S ou t- 

We also find that some groups change their location on the ve- 
locity plane in distant regions. The yellow squares in Fig.|3]indicate 
the cylindrical velocities of the LS groups. For instance, Hercules 
(group 4) in Sr is at the same cylindrical velocity as locally (peak 
inside the square), but in S ou t it has moved to lower v$. We see that 
this transition is continuously traced when we consider samples in 
small steps advancing from S ou t to Sr q , and progressively to inner 
Galacocentric radius. For Sin, Hercules does not appear as a group 
separated from the central structures. The v$ of the Wolf elongation 
follows a similar behaviour with radius. We also observe that some 
groups move in vr. For instance, Sirius (group 2) shifts to the right 
of its corresponding square (the one at highest V) for S ou t. 

The arrows in the upper left of each plot show the median 



directions of the radial (vi os ) and transverse (v t ) velocity. For these 
regions, as the radial velocity errors are small (e Vlos ~ 2 kms -1 ), 
the errors affect mainly the transverse motion (e„ t ~ 17kms~ 1 ), 
distorting the groups in that direction. For the case of S ou t, the 
velocity errors or distance biases would never move Hercules in the 
observed way. By contrast, for Si n a distance underestimation could 
shift Hercu les closer to the cen tre of the distribution. However, as 
showed in Burn ett et al. | <H)T 1), distances for giants may have a 
slight bias but in the opposite direction. 

We have evaluated the effects of neglecting the extinction in 
de riving distances for st ars in these regions. By using the dust maps 
of ISchlegel et alj J1998I) , we find an upper limit on the extinction 
in the J-band of 0.2 mag (by considering extinction at infinity) 
for the three subsamples, which would introduce a distance error 
< 10%, which is small compared to the 28% uncertainty due to the 
distance derivation method itself. We can also test that our results 
do not depend on the distance method by using only red clump 
st ars. We select t hese by J — K colour and gravity using the criteria 
in Sieb ert et al. I J201 la) and re-derive the distance from their J or 
K magnitude assuming absolute magnitudes Mj = —0.87 and 
M K = -1.55 dZ witter et alj|2008l) . This gives us subsamples with 
30-40% of the initial number of stars. Considering a dispersion of 
20% in the absolute magnitude for this population, we find that 
the transverse velocity errors are reduced substantially, from e vt ~ 
17 km s~ 1 to e vt ~ 11 km s~ 1 . Nonetheless, with smaller velocity 
errors and alternative distances, our results do not change. 

We perform another test to establish if errors only could ex- 
plain the differences in the distant velocity distributions. We study 
if these velocity distributions could belong to the same parent dis- 
tribution as the LS, given the velocity errors. To do this, we in- 
corporate the observational errors of Sr , Si n and S ou t in the LS 
velocity distribution. In practice, we use the proper motions and 
radial velocities of the LS but the distances and sky positions of 
the distant subsamples to generate 5000 simulated velocity distri- 
butions after error convolution. Then we compare the mock distri- 
butions with the original one (at Sin, Sout, or Sr q ), by computing 
the distribution of the Kolmogorov-Smirnov (KS) 2d statistic. In- 
stead of using this statistic directly, as a reference we also compute 
the distribution of the KS statistic between the mock subsamples in 
a given location themselves. The area of overlap of the two KS dis- 
tributions (obtained for the mock-original and for the mock-mock 
comparisons) thus gives us the probability that the distant veloc- 
ity distributions belong to the same parent distribution as the LS, 
given the velocity errors. Notice that our mock samples may have 
too larg e errors because firs t, the distance errors may be overes- 
timated jBurnett et al. I l201ll) . and second, the LS velocities have 
already their own error. We find that for the three subsamples, this 
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probability is < 0.06%. Repeating the same test but introducing 
also a distance bias of ±10, 20, 50 % does not increase the proba- 
bility. Only for Sin we find a 5% probability of consistency with the 
LS, given the velocity errors, if the distances were overestimated by 
10 or 20%. But in this case, this bias would move Hercules far from 
the centre of the velocity distribution, contrary to what we see. 

Fig. [3] (right panel) shows the velocity distribution of the sub- 
sample Z57 (median height below the plane of Z — — 0.6kpc). 
Again, some of the LS groups can still be recognized and are 3-cr 
significant. For instance, Hercules is here group 2, the Wolf elon- 
gation is present and group 3 has the same velocity as 7Leo. For 
a sample in the same cylinder as Z57 but with heights in the range 
-800 < Z < -600 pc (median of Z = -0.7 kpc) Hercules and 
the equivalent to 7Leo are still detected but with 2-er confidence. 



5 SUMMARY AND DISCUSSION 

We have analysed, for the first time, the velocity distributions in 
regions that are ~ 1 kpc far from the Sun using RAVE data and a 
structure detection technique based on the wavelet transform. We 
have found that most of the main local kinematic groups (Hyades- 
Pleiades, Sirius, Hercules, Wolf elongation, 7Leo) can be still de- 
tected at ~ 1 kpc from the Sun in the direction of anti-rotation. 
Some of them also persist at 700 pc below the Galactic plane. 
Therefore, these groups are actually large scale features. Given also 
the spread in age and metallicity of some of these groups (A08), this 
supports a dynamical origin instead of being dispersed clusters. 

In the distant disc volume at the Solar circle the groups ap- 
pear approximately at the same cylindrical velocities as locally. In 
contrast, for the distant regions (in the direction of anti-rotation) at 
different Galactocentric radii (±600 pc) the groups appear shifted 
in the velocity plane. These are the first significant changes in the 
kinematic structures ever observed in distant disc velocity distribu- 
tions. Especially at inner radii the relative configuration between 
the groups changes considerably. These trends are consistent with 
the typical spatial scale of variation of the kinematic structures in 
models of e.g. spiral arm effects (All). 

In particular, we show that the Hercules group moves to larger 
(respectively, smaller) azimuthal velocity for regions at inner (re- 
spectively, outer) Galactocentric radius, while it appears at the same 
velocity at the Solar radius. This is the behaviour expected for this 
structure if caused by the Outer Lindblad resonance of the bar 
(fig. 4 in [Dehnen 2000), although it can be also explained base d 
on the epicycle theory for a resonant feature dOuillen et alj|201 lh . 
We also observe shifts in Galactocentric radial velocity, which has 
also been observed in several models (All, lOuillen et al.l 1201 H ; 
iMinchev et al]20ld) . 

We have also found a new significant group in the Solar neigh- 
bourhood and confirmed the presence of structures that were not 
significant in other studies and samples. Some of the observed fea- 
tures such as the Wolf elongation are generally only vaguely con- 
sidered in the modelling. These and especially the observed groups 
in distant regions may constitute the keys for the understanding of 
the complex kinematics of the MW disc. In particular, the magni- 
tude and direction of the velocity shifts detected here in distant re- 
gions may tell us about the specific mechanism that produces them 
(bar or spiral arm related). Quantitative comparison with the mod- 
els could allow us to constrain the properties of the spiral arms and 
the bar of the MW, breaking the degeneracy found when fitting the 
models onlv to the main local erouns. 
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